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ABSTRACT 

11 Current sheets are ubiquitous in the solar wind. They are a major source of 
the solar wind MHD turbulence intermittency. They may result from non-linear 
interactions of the solar wind MHD turbulence or are the boundaries of flux tubes 
that originate from the solar surface. Some current sheets appear in pairs and are 
the boundaries of transient structures such as magnetic holes and reconnection 
exhausts, or the edges of pulsed Alfven waves. For an individual current sheet, 
discerning whether it is a flux tube boundary or due to non-linear interactions, 
or the boundary of a transient structure is difficult. In this work, using data 
from the Wind spacecraft, we identify two three-current-sheet events. Detailed 
examination of these two events suggest that they are best explained by the flux 
tube crossing scenario. Our study provides a convincing evidence supporting the 
scenario that the solar wind consists of flux tubes where distinct plasmas reside. 

12 1. Introduction 

13 The solar wind provides a natural environment to study MHD turbulence in a 

14 collisionless plasma. Over the past few decades the launches of spacecraft such as Voyager, 

15 Helios, Ulysses, Wind, and ACE have made available a significant amount of data for 
i6 analyzing the solar wind MHD turbulence. 



Kolmogorovl ( 1194 ll ) 



17 The first theory of hydrodynamic turbulence, suggested by 

is known as the K41 theory, predicted a magnetic field power-law spectrum ~ k~ 5 ^ 3 . This 

19 —5/3 exponent arises from the nonlinear interactions of the homogeneous hydrodynamic 

20 turbulence in which energy is cascaded from large scales to small scales. For incompressible 
2i MHD turbulence where the cascading process is mediated by counter-propagating Alfven 
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22 wave packets, the Iroshnikov-Kraichnan (IK) the ory (llroshnikov 



23 some recent theories of strong MHD turbulence (IBoldyrev 

24 predict a power spectrum ~ k~ 3 ^ 2 . 



2006 



.964; Kraichnan 1965) and 



Boldyrev and Perez 



2009| ) 



25 One important concept in turbulence is intermittency. iRuzmaikin et al.l (119951 ) 

26 suggested that intermittent structures can affect the solar wind MHD turbulence power 

27 spectrum. They pointed out that the effect of intermittency in the solar wind MHD 



2B turbulence is to reduce the exponent of the power law spectrum. 



Ruzmaikin et al 



f l!995h 



29 further suggested that intermittency is "in the form of ropes, sheets or mo re compl i cated 



30 fractal forms." Recently, in studying current sheets (CS) in the solar wind, 



Li et al 



(j201lh 



3i found that the power spectrum of solar wind magnetic field behaves as K41 in periods that 



32 have abundant numbers of current s 

33 current-sheet free (see also 

34 intermittency, the study of 



r eets a nd behaves as IK in periods that are almost 



Borovskv 



Li et al. 



20 



0)). Since these current sh e ets ar e a source of 



(120 111 ) supports 



Ruzmaikin et al. 



(I1995h 



35 A current sheet is a 2D structure across which the magnetic field direction changes 

36 abruptly. Current sheets can be of large scales. For example, the heliospheric current sheet 

37 and current sheets found in CME-driven shocks are all large scale current sheets. These are 

38 not the subjects of this study. Here we consider current sheets that are of small scales. 

39 Some current sheets occur in pairs. These can be tangential disco ntinuities (TDs), 



often forming the two boundaries of a magnetic hole (see the review of 



Tsurutani et al 



4i (j201ll ) and references therein), or rotational discontinuities (RD s) which are the bo undaries 



of an exhaust from a reconnection site (see (jGosling et al. 



20051 ) and the review of 



Gosling 



43 (120 111)1. Com paring to magnetic holes, reconnection exhausts can be of larger scales. 



44 iGoslingi (120071 ) has found that the typical width of a reconnection exhaust is ~ 10 4 km and 

45 some reconnection exhausts can be as wide as 10 5 km. Consequently, these boundaries may 

46 be practically identified as "single-current-sheet" event. 
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Most current sheets do not occur in pairs. These c urrent sheets can be generated 



throu gh non-linear intera ctions in the M 



2004J). Using ACE data, 



Vasquez et al 



TD t urbulence (IZhou et al. 



2004 



Chang et al 



(120071 ) examined magnetic field discontinuities 



which can have very small spread angles for Bartels rotation 2286 (day 7 to 33 in 2001). 
They found that the statistical properties of these discontinuities form a single population 
and they are consistent with turbulence generated in-situ. By examining the probability 
density functions (PDF) of the magnetic field components from a ID spectral code, 



Greco et al. 



PDF. Moreover, 



(I2008[) showed that cur rent sheets often occur at the super Gaussian tail of the 



Greco et al. 



(120091 ) found that, at the inertial scale, in which the energy 



56 cascading rate is independent of the scale, the PDF of waiting times ( WTs) between MHD 



discon tinuities that are identified in the solar wind using the method of 



Tsurutani fc Smith 



( 119791 ) and those from MHD simulations are very similar, suggesting that these structures 
can be explained as a natural result of the non-linear interaction of the solar wind MHD 
turbulence. 



Other opinions exist. In an earlier work, 



Bruno et al. 



(120011 ) studied current sheets in 



62 the solar wind by analyzing Helios 2 data using the min imum varian c e met hod to show how 



Bruno et al. 



(2001) were the first 



Borovskvl ( 20081 



63 the magnetic field changed over selected time periods. 

64 to suggest that these structures may be boundaries between flux tubes. 
es analyzed an extended time period of magnetic field from the ACE spacecraft and examined 
ee the distribution of the spread angle across the current sheets. He showed that the angle 

67 distribution has two populations and suggested that the second population, dominating at 

es large angles, could be "magnetic walls" and originate from the surface of the Sun. A solar 

69 wind that c onsists of many flu x tub es can be viewed as a structured solar wind. In both 



the work of 



Bruno et al. 



(120011 ) and iBorovskyj (120081 ) , the solar wind is envisioned to be full 



71 of structures, the flux tubes. Observed at a spacecraft, these structures are convected out 

72 with the solar wind. A similar scenario where structures convecting out from the Sun has 
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74 energy a r by 



Bruno et al 



73 been propose d by iTu and Marschl (Il99ll ). Analysis on the cross helicity ov and r esidual 



(120071 ) supported the proposal of 



Tu and Marschl (119911 ). 



Li 



(120081 ) devel oped a procedu re to 



Li et al. 



( 120081 ) examined 



75 Regardless the origin of a current sheet, 

76 systematically identify these structures. Using this procedure, 

77 current sheets in the solar wind and in the Earth's magnetotail using Cluster magnetic 

78 field data and concluded that current sheets are more abundant in the solar wind. Later, 



Miao et al 



( 1201 ll ) examined over 3 years' worth slow wind data using Ulysses observations 



and found there were 2 p opulations 



si sheets, in agreement with 



Borovskvl ( 120081 ) 



or th e distribution of the spread angle across current 



Perhaps a large fraction of current sheets identified in the solar wind are due to 
the non-linear interactions of the solar wind MHD turbulence, as shown in the work o 



Greco et al. 



(12008 



20091 ). However, a statistical study such as 



Greco et al. 



(12008 



2009|) 



as can not rule out the po ssibility that som e current sheets in the solar wind are boundaries 



se of flux tubes. Indeed, iBorovsky! (J2008I ) has used plasma data including proton density 

87 and temperature, Helium abundance, electron strahl strength, etc. to identify possible 

as plasma boundaries. Plasma data, however, is often of lower time resolution than magnetic 

89 field data. Furthermore, plasmas in different flux tubes may have similar properties except 

90 different velocities and magnetic field directions. Therefore, to unambiguously separate 

91 these two populations can be hard. Note, the occurrence rates of these two populations 

92 may have different radial dependence and/or different solar wind type dependence. 

93 In this work, as an effort to identify flux tubes in the solar wind, we present a case study 

94 of two "triple-current-sheet" event using data from spacecraft Wind. A triple-current-sheet 

95 event is where three current sheets occurred in a relatively short period of time. The reason 

96 that we want to search for a triple-current-sheet event is the following. In the flux-tube 

97 scenario, the solar wind plasmas reside in different flux tubes and the solar wind magnetic 
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98 field and plasma properties differ in these flux tubes. Since flux tubes are 3D structures, we 

99 expect the boundary between two adjacent flux tubes to be curved and have small-scale 
ioo ripples. This is shown in the cartoons in Figure [TJ As these flux tubes are convected out 
ioi past a spacecraft, depending on the relative configuration of the spacecraft trajectory and 

102 these ripples, one expects to observe most often a single crossing as in Figured] (a), and 

103 sometimes a double crossing as in Figured] (b), and occasionally a triple crossing as in 

104 Figure [1] (c). These three different cases are referred to as "single-current-sheet" events, 

105 "double-current-sheet" events, and "triple-current-sheet" events in this study. 

we A triple-current-sheet event can be used to discriminate between the scenario where 

107 current sheets are generated in-situ and the scenario where current sheets originate from the 

los surface of the Sun. In the former case, one expects no correlations between these current 

109 sheets in the sense that plasmas before and after these current sheet crossings need not 

no show any relationships. In the latter case, however, the spacecraft traverses through two 

in distinct plasmas in the sequence of "I, II, I, II" , so the observed plasma properties do not 

ii2 vary arbitrarily. 



2. Data Selection and Analysis 



We use the 3 s plasma and mag netic field data from the 3DP (ILin et al.l (119951 )) and 



us magnetic field (ILepping et al 



(119951 )) experiments on the Wind spacecraft. The data 
lie period was from September to October of 1995, which was during the declining phase of the 
H7 solar cycle. It is ideal to select data in the solar minimum period due to lack of transient 



us structures such as CM 
n9 readers are referred to 



s. For t he d ata analysis metho d of current sheet identification, the 



Lil (12008T ) and 



Miao et al. 



(120 llh . 



In the following, we first present a single-current-sheet event and a double-current-sheet 
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Fig. 1. — Cartoons illustrating (a) a single-current-sheet event, (b) a double-current-sheet 
event, and (c) a triple-current-sheet event. Note that in the case of triple-current-sheet 
event, the spacecraft traverses through two distinct plasmas in the sequence of "I, II, I, II". 
Consequently, we expect to find the plasma properties to vary accordingly. Dashed line 
with arrow represents the relative trajectory of the Wind spacecraft passing through the flux 
tubes. 

121 event. We then present two triple-current-sheet events. 

122 Figure [2] is a single-current-sheet event which occurred on 1995-09-21. The current 

123 sheet in Figure [2] is located at 14:32 UT and is shown by the brown vertical line. Before 

124 the current sheet the magnetic field magnitude \B\ decreases and the proton density 

125 N p increases. In the scenario where a current sheet is the boundary of a flux tube, 

126 these changes across the current sheet occur because plasmas in different flux tubes have 

127 different properties (as shown in panel (a) of Figured]). However, one need not invoke the 

128 flux-tube-crossing scenario to explain Figure [2j It can be simply a tangential discontinuity 

129 (TD) or one side of a reconnection exhaust. Indeed, careful examination shows that there 

130 was another small current sheet on ~ 14 : 32 UT, when \B\ decreased and N p increased. 

131 Our selection procedure did not pick out this earlier current sheet. 

132 The change across the current sheets are Alfvenic. The angle between 5B and 5V across 

133 the current sheet is 7°. For the earlier current sheet (that did not get picked up by our 

134 procedure), the angle is 173°. Such parallel and anti-parallel Alfvenic changes are always 



135 associated with a reconnection exhaust (JGosling||201ll ). Furthermore, there was also a 

136 decrease of magnetic field and an increase of number density (but not temperature) between 

137 these two current sheets, providing another support for identification of a reconnection 

138 exh aust. Therefore, Figure El although identified as a single-current-sheet event using th e 



2005 



Lil fJ2008r) algorithm, is re ally part of a pair of a bifurcated current sheet (IGosling et al 



iGosling 



2007 



20111 ). 



Figure [3] is a double-current-sheet event. Two current sheets can be identified around 
11:11:48 UT and 11:12:26 UT. In the scenario of flux-tube-crossing, Figure El can be 
understood as the spacecraft briefly crosses the magnetic wall between two flux tubes and 
then returns to the original flux tube. The schematic of this event is shown in the second 
panel of Figure [TJ Note that the temperature decreases and the proton density increases 
between the two current sheets. As in Figure |2} although Figure [3] can be explained by the 
flux-tube-crossing scenario, it need not to be. The angles between 5B and 5V across the 
two current sheets are 175 and 167 degrees. Unlike the first event, this double-current-sheet 
is not associated with a reconnection exhaust. There was a slight but insignificant drop of 
the magnetic field magnitude, so it is unlikely to be a magnetic hole. The proton number 
density and proton temperature were also changed slightly at the two current sheets. These 
slight changes, together with the pulse-like changes of t he 3 magnetic field compon ents, 



153 suggest that this structure could be pulsed Alfven wave (IGosling et al 



2011 



Gosling et al. 



20121) . No te. 



J2OI2J) 



154 if this was a pulsed Alfven wave, then according to Figure 3(a) of 

155 and the fact that it has a duration of 46 seconds, it would be a long duration pulsed Alfven 

156 wave. 

157 Figure |4] and Figure [5] show two triple-current-sheet events that occurred on 1995-10-2 

158 and 1995-10-13, respectively. Consider first the event shown in Figured! Throughout the 

159 event both B y and V y did not change much. B x underwent a sharp change at 22:13:05 UT, 
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leo the first current sheet, after which it only changed slightly, until 22:25:00 UT when another 

lei sharp change occurred and B x returned to values similar to those before 22:13:00 UT. At 

162 22:25:00 UT at the third current sheet another sharp change in Bx occurred. After the 

163 crossing, B x at and after 22:25:32 UT returned to a value comparable to B x at 22:24:30 

164 UT, just before crossing the 2nd current sheet. From the third panel, we can see that the 

165 V x changes at the same times as B x . 

lee Similar behavior also occurred to B z and V z . Before crossing the first current sheet 

167 at 22:13:05 UT, B z (V z ) was almost a constant. After the crossing, Bz increased and V z 

lea decreased. B z also became slightly more turbulent. At the 2nd crossing at 22:25:00 UT, B z 

169 and V z changed back to almost the same value as before 22:13:00 UT. Then both B z and 

170 V z underwent another sudden change at the third current sheet crossing at 22:25:30 UT. 

171 After the third crossing, B z returned to a value similar to that before the 2nd crossing at 
m 22:25:00 UT. 

173 The magnitude of the magnetic field \B\ (the 2-nd panel), the proton number density 

174 N p and the proton temperature T p (the 4-th panel) did not vary much throughout the 

175 event. Before the crossing of the second current sheet, around 22:24:30 UT, \B\ increased 

176 and Np and T p decreased. To better illustrate how the magnetic field direction evolves in 

177 this event, we have constructed an animation of the evolution of the unit magnetic field B. 

178 Two facts worth to note: 1) various plasma properties, including N p , T p and the 3 

179 components of B and V in the short period between 22:25:00 UT and 22:25:30 UT are very 
lao similar to those prior to 22:13:00 UT, suggesting that these are the same solar wind plasma. 
i8i This can be clearly seen in the online animation. 2) similarly, the solar wind before and 

182 after the short period are likely the same and it is different from that in 1). 

183 One may attempt to explain this triple-current-sheet event as the spacecraft crossing 

184 three uncorrelated individual current sheets that are generated by independent non-linear 
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185 interactions of the solar wind MHD turbulence. However, since independent current sheets 

186 have no correlations, the chance of the solar wind returning back to its original state after 

187 two independent current sheet crossings would be minute. Alternatively, one may argue 
las that the plasma between 22:13:05 UT and 22:25:00 UT represented a rather long-lived 
189 transient structure, and interpret the first two current sheets as the boundaries of this 
wo structure. In such a case, one has to explain why after the third current sheet crossing, 

191 both the magnetic field and the plasma return to values the same as inside the transient 

192 structure. 

193 Another triple-current-sheet event is the 1995-10-13 event, which is shown in Figure [5J 

194 Unlike the 1995-10-02 event, the 3-components of B and V, in particular, B x and V x , 

195 suffered some additional changes at and around the three current sheets, making the 
we 1995-10-13 event less convincing than the 1995-10-02 event. 

197 The first current sheet located at 19:32:28 UT. Both B x and B z showed a sudden jump 

198 across the first current sheet; V x and V z did not show significant changes. B y and V y also 

199 did not vary across the first current sheet. The current sheet is therefore non-Alfvenic. 

200 After crossing the first current sheet, B z was almost a constant for the next ~ 7 minutes 

201 until 19:39:40 UT, where the second current sheet was encountered. It increased across the 

202 second current sheet to a value similar to those prior to the crossing of the first current 

203 sheet. Comparing to B z , B x was nearly constant after crossing the 1st current sheet for ~ 3 

204 minutes and then gradually increased until 19:39:00 UT, after which it increases noticeably 

205 before the second current sheet. Across the second current sheet, it dropped to a value 

206 similar to those prior to the crossing of the first current sheet. The third current sheet 

207 occurred at 19:40:15 UT. Across the third current sheet, there was a significant change 

208 of B z and a small change of B x . The two black horizontal dashed lines indicate that B x 

209 (B z ) before the first current sheet was similar to B x (B z ) between the second and the third 



- 11 - 

210 current sheets. The two magenta horizontal dashed lines indicate that B x (B z ) between the 

211 first and the second current sheets was similar to B x (B z ) after the third current sheet. 

212 Note that the change of B x at the third current sheet was smaller than that at the second 

213 current sheet. After the third current sheet, B x kept increasing, until 19:40:30 UT. The 

214 value of B x after 19:40:30 UT is similar to those before 19:39:00 UT. As in the 1995-10-2 

215 event, we also constructed an animation of the evolution of the unit magnetic field B for 

216 this event. 

217 For the 1995-10-02 event, the angles between SB and 8V across the three current sheets 
2ib are 179°, 176°, and 174°, respectively. For the 1995-10-13 event, the angles between 5B and 

219 5V are 155°, 124° and 173°, respectively. While the three current sheets in the 1995-10-02 

220 event are highly Alfvenic, those in the 1005-10-13 event are not. 

22i 3. Discussion and Summary 



Current sheets are ubiquitous in the solar wind. They c an be generated in-situ t hrough 



223 non-linear interactions of the solar wind MHD turbulence (Greco et al. 



224 represent the boundaries of fl ux tubes that originated at the Sun (iBruno et al 



Borovsky 



2008 



Li et al 



226 reconnection exhausts (((Gosling 



2008 



2008 



2009), or 



2001 



Appe aring in pairs, they could also be the boundaries of 



201lh ) 



227 An intriguing question one may ask is: for any particular current sheet, can we identify 

228 how it originated? 

229 If current sheets that are generated in-situ and those that are convected out from 

230 the Sun have similar properties (such as the spread angles, the current sheet width, etc), 

231 then discriminating these two scenarios can be hard. However, as shown in the rightmost 

232 cartoon in Figure HJ the presence of triple-current-sheet event provides a strong support 
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233 to the flux-tube scenario. This is because in the flux-tube scenario the plasma and field 

234 changes across the three current sheets are intimately correlated: as the spacecraft crosses 

235 the three current sheets, the plasma before the first crossing and that between the 2nd and 

236 the 3rd crossing are the same; the plasma between the first and the second crossings and 

237 that after the third are the same. This is in stark contrast to the scenario where the current 

238 sheets are generated in-situ. In the latter scenario, the plasma changes at the three current 

239 sheets in a triple-current-sheet event need not match. 

240 Note that the identification of a triple-current-sheet event does not tell us how many 
24i single-current-sheet events are due to flux- tube crossing . As discussed earlier, since a 



242 reconnection exhaust can be of large scale ([Gosling 



20071 ). so me sing! 



e-current-sheet events 



Goslingl ( 120101 ) identified an 



243 we identify can be the boundaries of reconnection exhausts. 

244 occurrence rate of 40-80 reconnection events per month in solar minimum. In our study, we 

245 only consider current sheets which are abrupt (width <10 seconds) and whose spread angles 

246 are larger than 45°, we find about 350 "single-current-sheet" events per month. Assuming 

247 2 * 60 = 120 are boundaries of reconnection exhausts, then the rest are presumably either 

248 generated in-situ or are the boundaries of flux tubes. Assuming 80% (50%) of the rest 

249 are generated in-situ, then one gets about 45 (115) single- current- sheet events that are 

250 flux-tube crossings per month. 

251 If current sheets are boundaries of flux tubes that have a solar origin, e.g. super 

252 granules, then one may expect to find some statistical correlations between in-situ 

253 observ ations of current sheets and solar observation of super granules. Indeed, 



Bruno et al. 



254 ( 120011 ) have suggested that the sizes of the flux tubes, when tracing back to the solar 



255 surface, may corr elate with the size of photospheric magnetic networks. In the work of 



Miao et al 



( 1201 ll ). using Ulysses observation, the distribution of the waiting time statistics 



257 of the current sheets were obtained. Assuming these flux tubes do not split or merge during 
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25B their propagation to 1 AU, then one may expect such waiting time statistics resembles the 

259 distribution of the magnetic network sizes. Examining the waiting time statistics of current 

260 sheet, and in particular, its dependence on heliocentric distance, and its correlation with 
26i supergranule size will be reported in future work. 

262 To conclude, we have examined 2-month's worth solar wind data from the Wind 

263 spacecraft and identified two triple-current-sheet events. The sequence of the observed 

264 magnetic field and plasma data in these two events are in agreement with the scenario 

265 where current sheets are flux tube boundaries, as depicted in Figure [Q Unambiguous 

266 identification of flux tubes in the solar wind is important because these structures present 

267 an additional source of solar wind MHD turbulenc e intermittency. Th ey can affect the 



268 power spectrum of the solar wind MHD turbulence (ILi et al. 



269 the transport of energetic particles in the solar wind ( 



2011 



Qin and Li 



2008|). 



20121) as well as affecting 
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Fig. 2. — A single current sheet event occurred in 1995-09-21. Shown from top to bottom 
are the 3 components of the vector magnetic field in the Geocentric Solar Ecliptic (GSE) 
coordinate system, the magnitude of magnetic field, the 3 components of the vector proton 
velocity in the Geocentric Solar Ecliptic (GSE) coordinate system, and the solar wind proton 
number density and proton temperature, respectively. The brown vertical line marks the 
location of the current sheet. 
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Fig. 3. — A double-current-sheet event that occurred in 1995-09-26. Shown from top to 
bottom are the 3 components of the vector magnetic field in the Geocentric Solar Ecliptic 
(GSE) coordinate system, the magnitude of magnetic field, the 3 components of the vector 
proton velocity in the Geocentric Solar Ecliptic (GSE) coordinate system, and the solar wind 
proton number density and proton temperature, respectively. The two brown vertical lines 
mark the location of the current sheet. 
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Fig. 4. — A triple-current-sheet event occurred on 1995-10-02. Shown from top to bottom 
are the 3 components of the vector magnetic field in the Geocentric Solar Ecliptic (GSE) 
coordinate system, the magnitude of magnetic field, the 3 components of the vector proton 
velocity in the Geocentric Solar Ecliptic (GSE) coordinate system, and the solar wind proton 
number density and proton temperature, respectively. The three vertical lines mark the 
location of the current sheet. Also see the online animation of the evolution of the unit 
magnetic field B in this event. 
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Fig. 5. — Same as Figure HI but for the 1995-10-13 event. Shown from top to bottom are the 
3 components of the vector magnetic field in the Geocentric Solar Ecliptic (GSE) coordinate 
system, the magnitude of magnetic field, the 3 components of the vector proton velocity in 
the Geocentric Solar Ecliptic (GSE) coordinate system, and the solar wind proton number 
density and proton temperature, respectively. The three vertical lines mark the location of 
the current sheet. Also see the online animation of the evolution of the unit magnetic field 
B in this event. 
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